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ABSTRACT 


This work investigates the effects of digital processing 
in radiometers. It deals mainly with two digital versions 
wamamcoval power radiometer. | lhe first consists of RF, 
Mixer and IF sections followed by an analog to digital 
Ponverver., Sealimturther processing is done in a digital 
computer. The second version consists of RF, Mixer and IF 
sections followed by a square law detector, RC filter and 
analog to digital converter. From this point on the pro- 
cessing is done by a digital computer. A figure of merit 
is defined based on the performance of an analog total power 
radiometer. Exact results are obtained for the figure of 
merit of the first digital version. For the second, an 
approximate solution is obtained. The effects of saturation 
and finite step size of the quantizer were taken into con- 
sideration for the above results. The performance of digital 
balanced-Dicke and noise-adding radiometers is investigated 
using. the above results. The effects of digital filtering 


on the performance of a radiometer is considered. 
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ie INTRODUCTION 


A radiometer is basically a highly sensitive and stable 
noise receiver. The principle of operation lies in the fact 
that any object at temperatures above absolute zero radiates 
energy in the form of electromagnetic waves. Until some 
years ago, radiometers were used mainly as radio telescope 
receivers and constituted the main tool of Radio Astronomy 
[Ref. 1]. Lately, their use has expanded to many other 
areas. The measurement of atmospheric temperature [Refs. 2 
and 3], detection of air turbulence in the Troposphere 
[Ref. 4], airborne mapping [Ref. 5], the measurement of 
absolute radiation from a projectile flow field [Ref. 6], 
and passive detection [Ref. 7] can be cited as examples. 

The type of processing involved in a radiometer is 
especially fitted to digital metnods. A number of insti- 
tutions have integrated the radiometer with the digital 
computer in their “applications. The purpose of this work 
Ho vominvyesticate the effects of this integration, and to 
provide guidelines to the design of digital radiometers. 

The main body of the investigation will deal with two 
digital versions of a total power radiometer (TPR) shown in 
Figure 1. The first one (Figure 2) consists of RF, Mixer, 
and IF sections, followed by an Analog to Digital Converter 
(ADC) section. All further processing is done in a digital 


computer. The second one (Figure 3) consists of RF, Mixer, 





and IF sections followed by a square law detector, a low 
pass filter, and an ADC. From that point on, the processing 
is done by a digital computer. 

ic iid eeiwene PrenCemomesthe two cases 15 the location of 
the ADC. The second version (Figure 3) is presently being 
used and its main advantages over the first version (Figure 
2) is that it requires a much lower sampling rate. 

The results obtained for a total power digital radiometer 
Gan be used to determine performance of other types of 
radiometers. Performance of digital Dicke and noise-adding 


radiometers will be considered. 
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Pier ehtORMANCE HEGURE OheA RADIOMETER. 
DEGRADATION FACTOR 
By definition, the sensitivity, or minimum detectable 
temperature, of a radiometer is taken as that change in 
operating temperature T that causes a change in the expected 
value of the output equal to one standard deviation of the 
output [Ref. 1]. Refer to Figure 4. It shows the expected 
value and the standard deviation of the output voltage as 


functions of temperature for a typical radiometer. 
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FIGURE 4 
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From Figure 4, it follows that (to first order) 





sail = ee coe, Ae (1) 
op 


where E[-] = Expected value of [+], and T, 


D PomOpe ra vane 
Temperature. But by definition 


AE[I] =o = cert) 








ie Ca | T (2) 
Op op 
Therefore 
ae: 
AT = GETIT ja ) 
aT Op 
Me periormance’ figure can "be defined as, 
eT 
AT ar 
—S— O (4) 
Top TET 





JC, 
Op 

It can be shown [Ref. 1] that for an analog total power 
radiometer (Figure 1), using the fact that the inverse of 


the integration time is much smaller than the bandwidth of 
Piece bes Mixermmh filter, thar 


. - 2 
EL.) = tkT Sf JH (2) | af Cp) 


= ul 
Oyo = atk°T i gl TESS), eh aaseli (6) 





where 


Te— intvegration time 
Ho (f) = Band-pass filter (RF,Mix,IF) frequency 
response | 
k = Boltzmann's constant 


Bupscticuting into (4) 


: we SS (7) 


1. 
Top (B,t)* 


where 


( os [Ho (£)|° ary? 
0 


SS SS (8) 
4 
POE)” Ga 
7 | 


The performance figure of the analog TPR, (7), will be 
used as a basis for evaluating the performance of a digital 


TPR. Hence, the Degradation Factor can be defined as 


Wee rer lOrmoMeemuLsUromO lm Diettal nadLometver 


Performance figure of Analog Radiometer 


Intuitively, F should be equal or Proscar than one since 
as the sampling rate of the ADC goes to infinity and the 
step size of the quantizer goes to zero, both digital 
radiometers are reduced to an analog TPR. 

some comments are appropriate at this point. First, 


(5) and (6) show that for analog radiometers E[I] and oF 
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apeweangar LuUnceLoOnSs Of Pemperature "and? go through the 


origin. Hence equation (4) reduces to 


‘aigh 
a 
Op 


fe 

= sec (9) 

The addition of a quantizer, due to saturation, makes 
ELI] and or nonlinear with respect to temperature, and the 
performance figure is thus temperature dependent. There- 
fore (4) must be used instead of the better-known (7). 
Secondly, the parameter a in equation (4) is not a 
SOnvenleny cholce, Since i orgaggiven CE Gitferent gain 
SCLC Seem ENCenr ams Mixer, IF filter will change the 
performance of the radiometer. Referring to Figures 2 and 


3 it is easily shown that 
AE =e is! ve] (Figure 2) (10) 
(Figure 3) Geb 


where @ and 8 are proportionality constants. Using the 

parameters Jy 2 or To as appropriate, instead of T, norma- 
© 

mizaes the problem of gatingsetting out of the performance 


figure. 


US 





inehaecor ene rPediometer of Fiigure 2 
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Oy 2 
ACh Cc 
iT a | det (12) 
Op 5 
doy ~ 19x 2 ot 
C c op 
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AT _ z 
hae) GEL (13) 
P do co. atT . 
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III. DEGRADATION FACTOR CALCULATION WITH NO QUANTIZER 


An ADC performs two basic operations on a signal. One 
is sampling, and the other is amplitude quantization. In 
this chapter the effects of sampling alone are considered. 
Hence it is assumed that each sample can be represented 


without any error. 


A. IF SAMPLING 

Refer to Figure 2. Assume that at the antenna white 
gaussian noise is present with double-sided power spectral 
density No/2 watts/hertz. Then X(t) is a Narrow Band 
Gaussian Process with known power spectral density, and 


can be represented as 
X(t) = X,(t)cosenf t _ X(t )sinemf t uy 


where X(t) and X(t) are baseband Gaussian processes and 
Po is the IF center frequency of the radiometer. It can 
be shown [Ref. 8] that if X(t), X(t) and X,(t) are to be 


wide sense stationary, then 


Ry (u) = Ry (u) (15) 
Cc S ; 


Ry (yu) = ee Ei: + R 


x X (u)sinemf ou (26 } 
se 


ala 





where Ru) is the autocorrelation function of the process 
X(t) and Ry x (u) is the crosscorrelation function of the 
‘processes x(t) and X(t). 

Furthermore if H(t) is symmetrical about f (this As 


the only case considered in this work), then 
Ry y (uw) = 0 (17) 
cs 


and Ry (u) is given by 
C 


co 


Ry (u) = a S,(f)coslen(f-f )ulat (18) 


and Sy (f), the power spectral density of X(t), is the 
neurite, Transform of equation (18). 

Physically, for the case of interest (i.e. X(t), X, Gay) 
and X. (t) wide sense stationary, and Hy (f) symmetrical about 
myy4s 1 es is just twice the low pass version of S yf). 
ae to Figure 2, the low pass filters filter out the 


double-frequency terms coming out of the multipliers. Then 
IT=I, +1 CAD ) 

Ee 

x X, (kA) (20) 


since X (kA) = q, (kA) and 
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gt 2 , 
r X. (kA) (21) 
k=1 


since X (kA) = q, (kA). 


Hi 
He 
> [- 


Using the fact that for Gaussian random variables 


ELX, XXX, 1 = EC XX, JHK,X, JEL XX, JE(X,X, J+ECX, X, JELX5X51 





(a2) 
and expanding the double summation, it follows that 
4 
BX Af 
Or s = == es z (1-—)Ry (kA) (23) 
a k=1 Cc 
and 
E[I.] =o : (24) 
1 X 


Cc 


For the random variable I, the same results are obtained 
since i = Saget 


Then 


2 


X 
Cc 


E{I,J= 20 (25> 


and 


G23 = G6 GC + 2E(I,1,] = 2E[T, JECT, ] 
1 2 (26) 


iE) 





Expanding E{I,1,] and using the fact that R = 0, 


x x 6H) 
ons 
the last two terms of (26) cancel out. 


‘Using equations (25) and (26) in (12) it follows that 


~1 
2 @ al —_— Vetere 
(a = = Ir 2 (1 - ~) Ry (kA) (27) 
op no, k=1 e 
Cc 
or 
\ n-l 
Ae ee Ai OA kK 2 
et eee 4B Cee —) ROCA) (28) 
op TOy k=] e 
Cc 
where 
TS say (29) 


The degradation factor squared is 


CBA n=l 
Zi E Lae 
Fos BoA + —S- OE (1 - SE)Ryo (kA) (30) 


Ry (eis ff function of H,(f), therefore (30) will be used 
one | 
with four different frequency characteristics. 


1. Ideal Bandpass Filter 


Refer to Figure 5 


N, B sinmBu 


— (31) 


soe (CU) 
KX 
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=) 
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N jo 


MAM EUNSID, 5 
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and 


xX = Ne B oe) 


SUDSPuLUen es mest.) and (32j8into (30) 


rl 


Petes 2 


as sin’ TkAB 
E E a 


yen(i“ kt 


ry 
ti 


Oye) 


But 


Bo = BSW (34) 


where B, is the equivalent bandwidth from (8) and W is the 
half-power bandwidth. 


Then (33) can be rewritten as 


rae IL 2 2 jee 2 
= WA +°-2WA £ san os at > ksin 1 
k=l (1kWA) k=l (a7kWA) 


ne 


35) 

Since tT>>A, it can be shown that the third term of (35) 
is much smaller than the second and thus can be ignored. 
Letting A approach zero (sampling infinitely fast), (35) 
becomes 


15 O 
Re = OW ff Sin“mWt 4; (36) 


0 Ort )< 


and then since tT is large, the integral's upper limit can 
be replaced by infinity, and the integral reduces to one, 


giving themcorrect result. 
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Actual lyeeby Including the third term as A¥goes to zero, 
Paw iepeca slaty Lossschameone .= The reason for this’ts 
that in calculating the performance figure of an analog TPR 
iar, the muact. Ghat = << B was used to get (5) and (6). 
Hadmunls fact mov been used, aS "A approached zero, F would 
approach one exactly. 

Mewieiebe pletted as®a function, ef WAS Hence, for this 
case the Nyquist sampling rate (twice the highest frequency) 
corresponds to WA=1. 

The analysis done in this chapter is just a special case 
of the analysis that will be done in Chapter V, hence all 
results will be plotted in the figures of that chapter. The 
results obtained from (35) for an ideal bandpass filter are 
plotted in Figures 18 through 23. It is of interest to note 
that for W less than one (sampling faster than the Nyquist 
rate), F is exactly equal to one, an intuitively satisfying 
result. 


gree ingles luned Pilver 
For a Single Tuned Filter 


ae Jl, dl 
AS Oleg a Sea a oe 
1+4(—-2) 1+4(— =) 


where W is the half power bandwidth and W is much less than 


fos It follows that 


S,(f) = to IH (£) 4° (38) 
X ee O 


ZS 





and ‘ 


N 
oo (39) 
Xo 1+4(5)° 


From (39), using the Inverse Fourier Transform 





oe =n ul 
Exh = = a e (40) 
and 
AGEN AMEE KAD. —onkWA 
ae ae eee I (41) 
Op k=1 
Equation (41) can be put in closed form by letting 
Ae a 7eTWA (42) 
Using the fact that [Ref. 9] 
n-1 n 
=] L-A ? 
and 
n-1 n at sa | 
s kak = ee (44) 
k=] (1-A) 


and after some tedius but straightforward manipulations, it 


follows that 


| 3 | 
(Aa -)e = = (coth mWA) - ao eonccne TmWA)(1 - Q~emWT) (45) 
op 
and 
2 
ee (ABR) e 


etBy, 
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It can be shown that the equivalent bandwidth of a 


single tuned fulter is [Ref. 1] 


Be = TW (47) 


where W is the half-power bandwidth. Furthermore, since 
t>>A, the second term is negligible with respect to the 


ier soe rence 


F° = TWA(coth TWA) (48) 


To compare with the analog case, let A go to zero. 
Then F would become one, which is the correct result. 
Figures 24 through 29 show the results obtained for this 
filter. 
3. Gaussian Filter 


For a Gaussian Filter 


By (am Oe Se uaone 


Pry 9 = teh ete LN (49) 





where v = 41n2 = 2.773, Wis the half-power bandwidth, and 
Wis much less than fo: | 


It follows that 


~v(7)* 
hale = Noe (50) 
and 
T 2 
~—(Wr) 
Ty 45 
Ry (u) = NW(—) *e : (51) 


2D) 





Replacing equation (51) into (30) and using the fact 


that for this gfilver 


eye rt 
B, = (a = 1.505 W (52) 


itetollows’ that 


J n-1 ™ -£(mkwA)* 
F“ = 1.505WA + 3.01WA = OSS) Be (53) 


\ 
Using the same argument as before (A<<tT), (53) becomes 


; n-1 -£(mkwA)* 
re © Seay — 2 Crieeenn ae (54) 
k=l 


To compare to the analog case let A approach zero. Then 
(54) becomes 
a -<(nkwWa)* 


7 2 away 7 A 
0 


(559) 
and since tT is large, the integral's upper limit can be 
Hep laccdmpy sintinsty panemene amteorale recucesete one, as 
it should. The numerical results are plotted in Figures 30 
BaArougn 35. 
4, Butterworth Filter 
The last eer to be considered is a Butterworth 


Filter of order two. For this case 


a(n? = 
LP Leas) 


Ih 
1+16( 77 ) 





where W is the half-power bandwidth, and W is much less than i 
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Then 


N 








S, (f) = ——— 7 (57) 
Ko 1 +16(=) 
and 
N.TWW --Iu mW 
Ry (u) = =5 e fu (2)% cos(H Te - r) (58) 
C 


For this filter, it can be shown that 


2 we AC es Ty (59) 





CA) 


Substituting (58) and (59) into (30) it follows, 
after some effort that 
5 n-l 


FX = 1.48WA + 2.96WA CE e 
k=] 


3 
-k1(2) Wary a sin(ken (2) 2WA 


2 n-l 3 
- 2.96 (AW)™ r ke kt (2) Ward + sin(kn(2)*WA] (60) 
ui k=1 | 


Again it can be argued that the third term is negligible, 


therefore 


n-l Be 
Fo = 1,U8wA + 2.96wa 2 e7X™62) WArd 4 cin ie (2)%wa)] (61) 
aa 


Comparing to the analog case by letting A approach 


zero, (61) becomes 


T 35 
Po = 2,96wA s e782) Wha 4 sin(m(2) @we) Jat (62) 
0 


oa 





Since t is large, the integral's upper limit can be replaced 
pomeniaihuyesand Lo reduces fas it should to one. 
The numerical results for this filter are plotted 


in Figures 36 through 41. 


B. RC FILTERING AND THEN SAMPLING 

Refer to Figure 3. This type of digital TPR was studied 
by Ohlson [Ref. 10] for the case where no quantizer is 
present. Only the results and its implications are presented 
Molex 


It can be shown that 


N 
E[Z] = ¢ = (63) 
and 
N 
Pe ee (2° oe (64) 
“ é oB_.t 
Bake 


where g is a system gain parameter, Ng Lsaas. defined Jon 


page 17, and ate is the time constant of the RC filter. 


Also 
ECCI] = E(Z] (65) 
and 
a. fy | oa t-l kA 2 


where Ry (h) is theautocorrelation function of the process Z(t). 
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Substituting (65) and (66) into (13), we have 


kA 
n-1 ae 
1 2 kA 
ey Oe eo Ee 
Op E RC k=1 
Defining 
a = e 4/tRe (68) 


and using equations (43) and (44) it can be shown that 
\ 








: -T/t A (Ay oe 
(l-e t/ RC ) ( esch ) 
Fe = = coth a - eee COORG ame URC (69) 
RC RC (t/t pq) 


Equation (69) is of the same form as equation (47). 
Note that in that case it was argued that the second term 
was negligible since TW was large (on the order of 10°), 
Here the same assumption can not be made since ban <8 
(i.e. the RC filter bandwidth is much smaller than the 
bandwidth of the RF, Mixer, IF filter). The first term of 
equation (69) reduces to one as A approaches zero. The 
Second term however, contributes negatively. Therefore F 
can be less than one. dhesreason  homeunise ics sim vie odds 
tion of an RC filter to the radiometer, which has the effect 
of lengthening the integration time t. As A approaches zero, 
a comparison is being made of two TPR, one with an RC filter 
after the square law detector and the other without one. 
The one with the RC filter will have greater sensitivity 
(smaller AT) due to a longer effective integration time. 


Therefore F can be less than one. 


ae 





Neve Verwmtt Coe eliects offsampling only are of interest 
(not the effect of the addition of the RC filter), the 
second term must be ai pemardea by assuming that tT/tao = © 
et See ipa nate ol tet hemitilteringsis done by thesintegrator 
ier Oolc oy NGhewnG ill Der. 

As explained in part A of this chapter, the results are 
presented in Chapter V. Refer to Figures 46 through 51. 

In summary, the effects of an ideal ADC (i.e. no 
saturation and step size equal to zero) was studied in this 
chapter. Its main importance is the fact that it provides 
a lower bound for the Degradation Factor and, that under 


CoriamimecOonditions..a correction can be added to these 


results to account for the finite step size of the quantizer. 
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ITV. TRANSFER CHARACTERISTICS OF A QUANTIZER 
WITH A ZERO MEAN GAUSSIAN INPUT PROCESS 
In this chapter the effects of a quantizer on a zero 
mean Gaussian input process are analyzed. Two cases are 
investigated. Refer to Figure 6. It is desired to find 
the autocorrelation function at point 2 given the auto- 
correlation function of X(t) and the knowledge that it is 
Gaussian and has zero mean. Several approaches have been 
explored for the case of a quantizer alone. The one sug- 
gested by Kellog [Ref.11] was used since it can be easily 
extended for the case of a quantizer followed by a square 
law detector, and also the equations that it leads to are 
in a form suitable for computer programming. Figure / 
shows two quantizers, one with even number of levels and 
Poemounecrcewith oad munbDer ole leve is. 


Let N be the number of levels of the quantizer. Then 


_ 2 Giggs NZ) 
SS yy tte 
and 
» elie = (NFS) /2 | 


oul 





X(t q(t) 





SQUARE 
LAW 


Dew BCTU 






x Ct) hace) 


FIGURE 6 


A. QUANTIZER ALONE 


The autocorrelation function of q(t) is 


Ry ty) = ELQ(x(t) JQCx(tt+u]} 


Let 
X(t) = 0, 
and : 
2 GAS) = a 
Then 
Ru) = FF Q(t, )@(e,)F (Cy sty5u)dz, az, 
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where 


2 2 
_(0,-2P;o5,5,+o, ) 
D D 
’ (1-p, ) 20, 
£(0, 3553) lS eS HS) 
Oy On (1-py as 
and 
Ry (u) 
Px = RO) TT) 


Noting that the dependence of Ro upon u is only via Py > 


the notation Ri (Py) will be used for brevity. Then (75) 





becomes 
N N } X, ; X, 
R. (p#) = 2 qed fi Geen) Cla = da 
qk, a er) Mee ee lueee 
4=1 j=l hes Xy_y 
RS) 
Let 
BS = b 4/0 (79) 
B - b5/Sy (80) 
then (78) can be written as 
X X 2 _ Ae: 
+ _j ~(8) ~2py 8) B5*Bo ) 
-) aa p 
N N : 2(1-9,°) 
Roy) a? Sa} q;4, I f : dB, d8., (81) 
1 jar x, Xx 2) 4s 
elgily rell 2M l= Dee) 
ox oe 





but 
2 oan 
B,  - 2B,B5Py + Bo” = 
Therefore 


2 2 





2(1 - p,*) 
and (81)\ becomes 

N ON “Oy 1 
Ry(oy) = 5 pty ¢ oT 

j=1 j=l yen, 

oN 

Now let 
Bo ~ PyB)- 


[2 (1-py*) y 


2 2 2 


2 2 


“g2Za 


2(1 = py°) 


2 
B — e 2 


X,  [am(1-py) 7" 


ox 


it follows that the second integral can be written as 


x /O. —- +B 
2 {erfl aae ] - erf[ A i} 
L2Qi-py ) I. [2(1-p, )J 
where 
a 2 
erf (a) = ie PLERN 
Tr @) 


3D 


(82) 


(83) 


48,48, (84) 


(85) 


(86) 


(87) 








If a new function is defined as 








N q X./0., — PB XO Gn = pee 
W(oys8,) = 2 #{ erf[ 4+ ] -err | = 4 ** |] | (88) 
j=l [2(1-p,") J [2(1-p, )J 
then (84) can be written as 
x 
i: 2 
(py) - : ar —— g 
Ralpy) = =e oe f I Oeeip mec. Sew dB CGO) 
ae” (2n)? a=1 7 X54 a . 
ox 


Dividing both sides of (89) by a 


expressions for X, and or it is found that a new parameter 


and substituting the 


appears, namely Oy /L . This parameter together with N and 
Py will be the independent variables of (89). Note that with 
the introduction of Oy /L » the dynamic range of the quantizer, 


ieee LS mo longer) an independent factor. 








pyle uu s 
qd; = a,/L (90) 
Xi = X, /L | (91) 
and 
| Oy = 9y/L | (92) 
ie Lol Lows. cha 
R_(p.,) . #,' 
p gt 6) om 
Ge x il ai x 2 
a = To Se ae 1 J wip ysB.) e dg (93) 
2. (Ene 5. eee ot 
| a 


and 
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Naa: X!/o.. - p,B a Ao} oN: 
¥" (py 5B, ) = ¥ =r {eret ee See ee = erf fe ae 1 (94) 


[2(1-p,°) 7 [2(1-p,°) 7 


It can be shown, following a similar procedure, that for 


Py equal to one, (93) becomes 


R(1) 62 N  q! '/o. xX! ,/o! 
i= a1. = 5 y (49° feret =i 1 - eret 44 | (95) 
Oy x? i=1 © (2)? (2)2 


This gives us 04° numerically so we can get the normalized 
autocorrelation of q. Equations (93) through (95) allow the 
computation of the autocorrelation function of the output of 
@a quantizer of N steps and saturation level L given that the 
input is Gaussian, zero mean, and of known autocorrelation 
function. The procedure would be as follows. hast, 
calculate O's Xs and qs for the given quantizer and input 
process. Second, find a value of Py fOumasparc.icular 2 
(Equation (76)). Third, form the function v' (py ,B,) for 
the particular value of Py determined in the second step. 
Fourth, replace the function v"' (py 58,) into (93) and perform 
the integration. This procedure will map one point of the 
input autocorrelation function into one point of the output 
autocorrelation function. 

A program was written to solve (93) through (95) for the 
case of symmetrical quantizers (See Appendix under Transfer 
Characteristics of Odd level Quantizer, and Transfer Charac- 


teristics of Even level Quantizer) and the results are 


Bi) 








plotted in Figures 8 through 14 as a function of Py- The 
Rg lu) _ 


ordinate axis is aD Py - in other words, for a given 
Py » Pq its obtained. The advantage of plotting it this way 
is that the curves are independent of the input autocorrela- 
tion function (as long as the input process is Gaussian and 
zero mean). 

In general, the curves are a function of 0,/L (except for 


the special case of N=2). Note that for N=2 the result is 


\ 
the well-known relationship [Ref. 7] 


Ry (Py) 
5 * 
q 


= = arcsin Py (96) 


It is of interest to note the behavior of these curves as 
0,/L approaches extreme values. For a quantizer with odd 
mumber®of levels, as Oy /L eets very large, the quantizer will 
appear as a hard limiter, hence the curves will approach that 
of Figure 8. As 0y/L approaches zero the quantizer will 
appear as a three step quantizer (eventually it will become 
an open circuit) and the curves will show a behavior similar 
tO those Shown in Figure 9. For a quantizer with an even 
number of Novela the behavior of Peale as 0,/L becomes 
large is also as described above. When 0,/L approaches zero, 


however, the quantizer will appear again as a hard limiter, 


therefore the curves will approach that of Figure 8. 
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B. QUANTIZER FOLLOWED BY A SQUARE LAW DETECTOR 
The reason for analyzing this case will become apparent 


in Chapter V. Refer to Figure 6 and observe that. 


¥(t) = q°(t) (97) 
and 

Ot ae QL ACG) | (98) 

Ry(u) = ELY(t)¥(t+u) 1 (99) 


\ 


Upon substituting (97) and (98) into (99) it follows 
that 
Ry (uw) = E{Q°LX(t) JQ°LX(t+u) J} (100) 


Following a similar procedure as done for the previous case, 


(100) becomes 








Xs 
= 2 
' ; -B 
Ry ( py) 1 Ne 4585 oX 1 
= 2 we 6G (101) 
lee ee ae NT! © 1 
oy (Or 2a =1 ext = 
Ox 
and | 


N q! X!/ol - pyB Xi _,/oy, - 9,8 
Mpys8,)=5 2 (Herel +444) - er HA i} (102) 
jel oo E2(1-p,")] [2(1-p,") J 


For Py=l it can be shown that 


FLY} _1 4 qi y Xi Xi 
Ts > (—-) [Lerf ( ——— ) - erf (—— ) ] (103) 
Of 2 yy Wey " ser 
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re qa» Oy, and py were GemineaParig(90), (91), (92) and (77) 
respectively. 

A computer program was written for this case when the 
quantizer is symmetrical (See Appendix under "Transfer 
Characteristics of Odd Quantizer Followed by a Square Law 
Detector" and "Transfer Characteristics of Even Quantizer 
Followed by a Square Law Detector"). The same comments as 
in the previous section apply for the plots of this case 
SmeCDomtuniat mowsthe ordinate akkiis ds 2a instead of Py 


since Y(t) is not zero mean. The results are shown on 


Figure 15 for N equal to five and eleven. 


StU, 
Y 
iGgecan Dey Shown that for the case of no quantizer Ry (0) 
is given by 
2 
Ryo) 3 sage 


Note how the curves of Figure 15 depart from the curve given 
by (104) due to the effects of the quantizer. 

In summary, in this chapter an algorithm has been devel- 
oped: suitable for ois programming in order to find the 
autocorrelation function of the output of a quantizer (or 
quantizer and a square law detector) given that the input 
process is Gaussian and zero mean, and with known auto- 
correlation function. What it basically does (for both cases) 
is to reduce a double integral (75) and (100) into a single 
integral by using the error function which is available in 


most computer libraries. 


a7 





The transfer curves for both cases, shown on Figures 8 
through 14 and 15 are independent of the input process 
providing it is Gaussian and has zero mean. These results 


will be used in Chapter V. 
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Vay pOBGRADATION FACTOR CALCULATION WITH 
BOTH SAMPLING AND QUANTIZATION 
iMimeciads chapverecme effects of slampling and quantization 
On the performance figure of the two digital radiometers 
under study are investigated. As in Chapter III, the two 


cases will be considered separately. 


A. IF SAMPLING 
Paeeoxacus RESULTS 
Refer to Figure 2. In Chapter IIT an expression was 
| Geveloped for the degradation factor for the case where there 
was no quantizer — see (30). Those results cannot be applied 
directly to this case since in order to obtain (30) it was 
assumed that q, (kA) and q. (kA) were samples of Gaussian 
processes — see (22). With the addition of a quantizer it 
is clear that this assumption no longer holds. Hence we 
must start again with a more general approach. This can be 
accomplished by concentrating first on the top channel of 
the radiometer of Eigurd 2. The processing that the ADC and 
computer do to the signal can be drawn sequentially as in 
Figure 16. It is argued that the order in which the signal 
qa, (t) is sampled and squared is immaterial as far as the 
random variable I, is concerned. Hence Figure 16 is a valid 
model of a channel of the radiometer of Figure eis 
' Referring to Figure 16 it follows that [Ref. 12] 
n-1l 


rE (1-S)(Ry(kA) - E°LYI] (105) 


Z 
Nel 
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Ry (k ) is the yen sample of the autocorrelation function of 
vat ) 


and 


2 


ELI, ] = Ely(t)] = o, (106) 
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Using (26) and the fact that for the case of Figure 


2 the RF, Mixer, IF filter is symmetrical about fo: then 


O a. 2 
pe we Sh (107) 
ell 
and 
“ ee 
2 = 2EL[T, J =? ae | (1:06) 


Cc 


Applying the definition of degradation factor it 


follows that 
2 


20 4BoA n=-1 

Ne kA 4 

=e A + ¥ (1 - —) [R,(kA) - 02] 

: aE of, k=1 Tk 4 

i a ee ee (109) 
2 : 
1 att ie "x, at Top 
do, 


Cc 
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Using the same argument as in Chapter III (A<<t) it 


can be shown that the term 
Tee RC AyEe Gael (110) 
gE v q 


can be disregarded due to its small value compared with the 
other terms in the numerator of (110). 


Then (109) becomes 





A+ > (Re(kA) = Gee | 
2 Sx a oe ae do 
Po = Guia), 
2 
doy Hee P 
Cc 


Looking at (111) term by term, it is noted that in 
the numerator, the first term can be obtained from (103) and 


(95). The second term can be obtained from the corresponding 
Ry, (kA) 
c 
malue of - and the t sfe haracteristics of Ba 
a O Rx, co) 2" ran r character @ 


quantizer followed by a square law detector (Figure 15). 


Actually, a-scale factor of oy? /oy will appear because of 
the way the curves were SR Eee The denominator can be 
obtained from (95). Note however that Equations (93), (95) 
and the transfer characteristics developed in Chapter IV have 
dy/L as a parameter and not a . This is due to the choiee 
of dy/L as a parameter since it was found in that chapter 


that the ratio of oy and L was of importance and not their 
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absolute magnitudes. Therefore it makes sense to have (111) 


be a function of Oy ig” eho Tan instead of Ty + Oy /G can 


be thought of as a beaiacta standard ney ie of X(t) and 
its value can be related to temperature in the same form as 
Equations (10) and (11). 

A plot of E[IJ/L®° versus oy /LE is shown in 
Figure 17 for different values of N. The effects of satura- 
tion of the quantizer can be observed as temperature (equiva- 
Iént to oS / i) is increased. 

The four filters previously considered for this case 
were used again and the results are shown in Figures 18 
through 41. The bottom curve for each figure represents the 
case investigated in Chapter HII where there is no quantizer. 
The programs used to evaluate (111) are included in the , 
Appendix. It should be pointed out that the programs used 
to calculate the transfer characteristics of Chapter IV did 
not provide a continuous function as it is shown on the plots 
of that chapter. Since in order to solve (111) the transfer 
curves must be known at every point, a least-mean-squared- 
error: polynomial Ral oy Ee done to the data obtained from the 
programs of Chapter IV. 

Looking at Figures 18 through 41 it is noted that in 
each case there is a reasonable optimum value of re 
where the degradation factor is lowest. Physically this can 
be explained by the fact that for small Teo , the step 
size degrades the performance of the system, hence increasing 


Jy /G will tend to minimize this effect. However a point is 
C 
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reached where saturation effects start dominating and the 
further increase of cae causes an increase of the 
degradation factor. For N>4 , raat optimum is around 
0.35 regardless of the sampling rate. This would imply that 
for a given + op: the gain setting of the REF, Mixer IF filter 
should be set at a value such that fe. = is equal to 0.35 
for optimum performance. Note that this optimality is for 
small source temperatures. For large signals, dynamic range 
is a eon ven WIE (Cloieyerr Sig \fdi 

é. An Approximation 

The evaluation of (111) for the exact value of the 
degradation factor requires extensive computer time. For a 
quantizer Sr LiVewbacs Oerreater they Computrarion, FCImMmes per 
Case is very long (on the order of hours). Therefore it is 
of interest to investigate an approximate solution for the 
degradation factor and to compare the results with the exact 
Pues already obtained. 

Figure 42 shows one of the channels of the radiometer 
on Ii gsuremc. Assume that the sampling rate is not too fast, 
the step size of the quantizer is small with respect to the 
standard Beacon of the process X(t) and no saturation 


takes place in the quantizer. Then 


q, (kA) = X (kA) + e (112) 


k 


where X, (kd) is the KB sample of X(t) and e, is a uniformly 
distributed random variable that accounts for the error 


introduced by the quantizer. The probability density 
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function of ey, is the following [Ref. 13] 


1 a a 
—. = ee < e < = 
f(e,) = + “ a (113) 
0 le, | 34 
k 2 
Then 
Yn Nn 
es eg tea] = £ [x2 (kA) + 2e X (kA) + ec} ~~ ~<C4) 
1 n = C n 5 Cc kc k 
k=1 k=1 
and 
male. Z i se eae 2D 
ELT, J —* )) ELX (kA) J + = ») EL2e,X, (kA) ] a 5 Ele; J (115) 
k=1 kK=1 k=1 
— ee ee ee 
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Assuming e, and X (kA) are uncorrelated, which is a 


good assumption if the step size is small compared with o 





Ko 
2 aé 3 | 
E{T, J = oa eo (116) 
Now 
= x 
etre} =e et (jA) + 2 E e,X,(JA) + : ae; 12 
n® j=l © ell Jeu 
Gale) 
n 
rE (ud) + 2 E e,X (kA) E ef J 
k=1 k=] k=1 

Expanding (117) and using (23) and the assumption 
that e, and e, are independent PrOreK ese) ,wwhichs is gcod if 
the sampling rate is slow, it follows that 
2 _ i, , 0 2 eye) i & 
= = A + a) z (n-k)Ry (kA) + = “3 Oy a ETO (118) 

1 n k=l e eG 

It was shown previously that for the case of Figure 

2) taiy Hate) was symmetrical about f_, then 
6 = 20° (119) 
I I 
al 
and 
Eft] = 2 ECT, J (120) 


al 





Using (119) and (120) it follows that 


Pek eo! Ff OS) ee wel Spa+ 2 Jens (121) 
Some’ rem n? By Jeo ay oe 
Oy k=l C Oy 0 
X 
Cc Cc Cc 
But the first two terms are just the degradation 


factor for the case of no quantizer, therefore (121) can be 


written as 


Ro = FC a. “A (122) 
. NQ 6 a E 360 a E 

X X 

Cc Cc 


where FQ is the degradation factor with no quantizer. The 

last two terms of (122) represent a correction, due to 

quantization, for the degradation factor with no quantizer. 
The approximate solution is very good if the assump- 


tions under which it was derived hold. For N=l17, oy = 0.3 
(a/oy = 0.42) and WA= 1.0 the error is less than one percent. 
For ties Of: A ale greater than 0.3 the approximation is not 

good due to saturation. For values much smaller, it again 


fails. due to the large value of a/O . Figure 43 shows how 
C 
the two results compare for N=7 and different values of 


Jy /L. In it, ‘the behavior of the two ‘solutions as oy /L is 
e ~ ae 
dnereased can be observed. For N much greater than 7, the 


approximation is very good when Oy /L is near 0.3. 
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B. RC FILTERING THEN SAMPLING 

Refer to Figure 3. Unfortunately there is no exact 
solution for this case since there is no way to get the 
second order probability density function of the process 
z(t). This would be needed to approach the problem in the 
same way as done for the radiometer of Figure 2. Two 
approximations will be discussed. 

1. First Approximation 


It was shown in Chapter III that 


ELZ] 


g(N/2) Gi2z3)) 


and 


2 2 ame 1 
gE (N/2)° sR (124) 


BE RC 


2 
7, 


Owing to the Central Limit Theorem it can be assumed 


that z(t) is a Gaussian process if B >> ieee Due Wo. the 


E'RC 
square law detector, it is obvious that z(t) can never be 
negative, hence the assumption is good only if the mean of 
z(t) is sufficiently large compared to its standard 
deviation so that the assumed Gaussian probability density 
function is’ negligible for negative arguments. 

It is also assumed that the quantizer has a bias 'b!' 
as shown on Figure 44. When b=0 this quantizer reduces to 
the one shown in Figure 7. It is argued that the order in 
which z(t) is sampled and quantized is immaterial as far as 


the output I is concerned. Hence Figure 45 is a valid model 


of the processing done by the ADC and computer. 
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FIGURE 45 
Let 
Elyza] = m Ci25y) 
then the probability density function of z(t) is 
_ (zem)* 
f(z) = e 207 (126) 
Z 2 35 
(2107) 


It can be shown, using a procedure similar to that 


introduced in Chapter IV that 











N , X47 X, .-m 
ELq(t)] = 2 2 Deng = )) SRST it ed (127) 
7 iM (2) On 2) "ay 
X,-m 
N On -g°/2 
Yor S Se GOR Gl ii P(pn58) e dg (128) 
a z (27)? 4=1 X,_,-m Pz 
oz, 
See = 0 B “oe! - 0 B 
N q On ja On Z 
Woy 5B) = Ze {erst | (129) 
i [2(1-95)]* [2(1-p5) 1" 
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‘But from Figure 44 


X, = 2L(4 - B)/(N-1) + bd (130) 


i 


and 


Gee CLT = (Nay) /2 17 (N-1) (131) 


Equation (130) can be written as 


X =X +b (132) 


where ,, is the value of X, when the quantizer is zero 
centered. 

Substituting (131) and (132) into Equations (127) 
ouch (129), the parameter b-m appears. Physically, this 
is the offset of the mean of the process z(t) with respect 
to the center of the quantizer. When no source is present 
the value of this offset will be determined by the quiescent 
conditions of the radiometer. When looking at a source, 
however, the value of the offset will change since the mean 
of z(t) changes linearly with temperature. 


Let 


eS iS fe CLS ey, 


where m1, is the expected value of z(t) at T. Then (127) 


a 
through (129) can be written as 
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i=l of (2) “On 
(134) 
; 55 t t 
es oe le a 
a) 
(2) *o5 
Xo nf 
\ i” 
7 Z 2 
R-(p7) Ne “B fe 
Ses 1, it logsB) © ap (135) 
O7 (27)? i=1 Z Xo.-1 
ope 
Z 
Xo ale One lee 
; N a o7' ae P78 eG te. a O7B 
b'(o758) = 2 {eret Ea ee SSG) le oe aE i} (136) 
j=l 203 [2(1-02)}* [2(1-08)}% 


where Oo mSephne value Oteipne standard deviation of Z(t) at 
O 


Top? and all primed variables imply that they have been 


normalized by L. 


For the case of P7 = 1 (135) becomes 


2 N a a fi 


Equations (134) through (137) can be solved in the 


same form as done in Chapter III. However, the addition of 
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1 
two new parameters, (f' and (2Btn0)?) makes the problem 
GQifficult to solve in a general way due to the number of 
variables involved. 


One way to get around this problem is to assume that 


1 
+6 
(2Bot pa) is large. Then 
N a 
ELq] = 2 qs if f,(z,m,o) dz (138) 
j{=1 x 
je] 
where 
_(z-m)* 
i 20° 
f,(z,m,o) = ——s—~ (139) 
7 oy 4s 
(210 ) 
From (123) and (124) 
8 mn 
“= O7 (2B, tao) (140) 
Therefore 
ap ON X,  df,(z,m,o) dz ; 
aims ” qs f SS (141) 
m $=] x dm 
jel 
but ; 
ar 
pia stm, af ao 
an om 30 am (142) 
Substituting into (141) and using the fact that 
of 
te ops 
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re LOMmowse® Chav 


dF Le N N X, BET 
a (2B tac) z q, [f,(X,_, sm,0)-f, (x, ie) a] se z Sie A 557 do (144) 
i=1 i=l 
i-1 
‘ 7 
1 
if ORS) is large, it means that the expected 


value of z(t) is much greater in value than its standard 
deviation. Therefore it can be assumed that a change in 
temperature affects the expected value of z(t) but the 
effects on On are negligible. In that case the second term 
of (144) can be disregarded. 


Using (66) it follows that 


a 2 2 
os: 6G n-1 R.(kA)_ = E-[q] 
=, Wig ere — a (145) 
oo Z k=1 On 


Substituting (144) and (145) into (13), it is easily 


shown that 





0 sane = Se ma Ry (ka) ~ Ea] 
7) s.. + Betac€ 2 
= eo? ie wig a ge 7 
| dELgq J 2 
do, 
2B, (A/ta.) Ss ae Ri (kA) - E Ca] 
ftpo el a cl 
iL 
dE{g | 2 
dg, 
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The parameter 2B.tp,, on the numerator of (146) will 
cancel out. The third term of (146) can be neglected only 


ff) T>>t (see comments in Chapter III. The addition of a 


RC 
quantizer does not invalidate them). 

It- must be pointed out that the parameter qe) *. 
even though it does not appear in (146), is an important one. 
It is a measure of how much smoothing the RC filter does. 

A large value makes quantization difficult. On the other 
hand, a small value invalidates the assumption that z(t) is 
a Gaussian process. Other consequences of this parameter 
will be discussed in Chapter VI. 

Equation (146) can be evaluated using (134) through 
(137) and (144). Basically it is done in the same form as 
for the case discussed in Section A of this chapter. Using 
(135) and (136) a set of transfer characteristics has to be 
computed similar to Figures 8 through 14 (that set of curves 
is the solution when f'=0). The major difficulty lies in 
the additional parameter f' which increases the dimensionality 
of the problem. 

The case when Teg (no offset) was analyzed using 
Figures 8 through 14 and (144). The results are shown in 
Figures 46 through 51. For any other value of f' (134) 
through (137) have to be evaluated. 

As an example the case of N=13, o7/L=0.2 and f'=0.5 
was solved. The results are shown in Figures 52 and 53. 

Note that they are very close to the zero offset case for a 


seven level quantizer with o7p/L=0.4. The reason is that 
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a/o 7 is the same for both cases. They are not exactly 
equal however, since in one case clipping occurs essentially 
for positive values of the signal, while in the other case 
clipping occurs symmetrically for both positive and negative 
values of signal. Figure 52 is a plot of the transfer char- 
acteristics for this example. The programs used are in the 
appendix (See "Transfer Characteristics of Quantizer With 
Offset," and "RC Filtering Then Sampling, Offset Case"). 

The case when N=2, f'=0 is of interest since results 
are obtained easily without extensive computer programming. 


Using (146) and (96) it can be shown that 





=. fe GR A n-1 (kA) oe n-1 R, (kA) ae 
= a eearcsin ~ % k arcsin 147 
Atpo tro 11 Ry (0) el R7 (0) 
RC 
The last term can be neglected when T>>t pa: As the 


Sampling rate increases, A goes to zero and Equation (147) 
becomes 


at 2 
Fo = f arecsine" ay (148) 
For large tT ;. the upper limit can be replaced by 


infinity and (148) integrates to 


2 


ie In 2 (149) 


ro} 


therefore 


YP oe ee (150) 
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The results are surprising since it means that a 
hard limiter degrades the system only by four percent. 
However, the dynamic range of a radiometer using a hard 
limiter would be very small hence this result is of little 
practical use. More on the problem of dynamic range will | 
be discussed in Chapter VI. 

2. Second Approximation 
Using the same arguments as in part A of this 


chapter, it can be shown that [Ref. 9] 


2 
Po = F ee ch. 


(151) 
NQ tre ou 





The results obtained from (151) are in close agree- 
ment with those obtained using the first approximation even 
for low values of N. It must be kept in mind that the 
second approximation does not take into account saturation 
effects of the quantizer. For the case when f'=0 it was 
Rouncechat if o//L < .6 WA>l , a/o, <2 and N> 5, 
the second approximation was within 1.5 percent of the first 
(For N=17 and the Rove conditions, the difference was less 
than 0.6 percent). These results indicate that a sizable 
amount of saturation can occur before (151) no longer holds. 
Figure 53 shows how the two approximation compare for the 
example done on page 91. Figure 54 shows a comparison for 
N=7, f'=0 and different values of o/L. The different 
behavior of the two approximations as o//L is increased can 


be observed. 
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To summarize, in this chapter an exact solution has 
been found for degradation factor of the radiometer of 
Pictucw ewe Or ule radiometer oferigure 93800 pmomcoineien 
exists, but an approximation can be made under the reasonably 
good assumption that z(t) is a Gaussian random process. For 
both cases the major difficulty in global analysis lies in 
the number of parameters involved. 

am approximate solution was developed which is in 
elose agreement with the first solution. Its major advantage 
is that it is easy to compute and as the number of steps 
finereases, i pets better. For practica® ADC it ishoudd be 


adequate. 


pa 


Midgets a 
' = 2 + 
° ae feet aess 
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=a 
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VI. DYNAMIC RANGE AND LINEARITY OF A DIGITAL RADIOMETER 


In this chapter the effects of the addition of an ADC 
on the linearity of theoutput of a radiometer is investigated. 
The main problem here is the saturation effect of the 


quantizer. 


A. IF SAMPLING 

Refer to Figure 2. The output voltage versus temperature 
characteristics for this case was investigated previously, 
and is shown in Figure 17. If the ADC was ideal (no satura- 
tion and a/o = 0) there would be one straight-line curve 
with slope equal to two. It can be seen that for small 
oy¢ /Lé (proportional to temperature) the curves tend to have 
a Slope of two. 

Figure 1/ provides all the information needed for the 
problem under investigation in this chapter. However, the 
way the curves were normalized makes them awkward to use. 


They can be changed starting with (108) and (95). 


Be 
: oO 
Xe Ne ig! 22 J Xe 
za = — 2 (=) [erf (== ) -erf (— L J Cle) 
ie L all Xo o,(2) 0¥(2) 
but 
Oo 
we 
Cc 2 
T° Weave (153) 
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Therefore (152) becomes 


2 
O N q,(N-1) 2 X, (N-1) (N~1) 
Bt a 4 y [= 5 ] fertt a -~ erf[ “t =~ 1| (154) 
a ye ah aoe 2(2) "oa 2(2) "a, /a 
e c 


Equation (154) is plotted in Figure 55. The curves asymptote 
to 


ELI] _ (N-1)° 
a nee GIUGIS » 





a 


This can be shown by using the fact thas as Oy/a 
gets large the quantizer appears as a hard limiter. 

Figure 55 shows that a substantial range of linear . 
output can be obtained with this scheme before saturation 
dominates. 

In order to use Figure 55 a criterion must be 
._ defined for how much departure from linearity is allowed 
for the temperature range expected (one percent, for example). 
Once that is defined, Figure 55 can be entered to determine 
the number of steps of the quantizer required. The value of 
the Away at Top will be determined mainly by how much 
degradation of performance is allowed (Chapter V). We thus 


have a trade-off of linear dynamic range versus minimum 


detectable temperature AT. 


INNS 





SOA ER 
NM | | 
EDGWARE 
ee S| 
NaN . 
BRRaNiii 
SA ELL 
RRA 
|| ja} | | AY 


FIGURE 55 
106 


K-SCALE=3.0Q0F+00 UNITS INCH. 
P-oURLE=5.00E200 UNITS INCH. 














B. RC FILTERING THEN SAMPLING 

Refer to Figure 3. A similar analysis can be made for 
this case. However the additional parameters (f' and 
eB tan) *) make a general representation difficult. 
Basically, (127) has to be rewritten in terms of Op/a 


instead of O7/L. It can be shown easily that (127) becomes 








1 
N aq! Xt = f+ (o! ~ol)s(2B_t.,) 7 
eae! fevet O ana: Bate 
a i=l 2soy (a)aora, 


X41 ~ £ + (of -09)s(2Bst a9) * 


erf[ . C5 i) 
(2) & ol 
vA 
where f', qi and X; are as previously defined, and 
7 
oy = o//a (15ife 
Glameece/a at” T (158) 
“s = 2/(N-1) 
Figure 56 is a plot of (156) for the case of f'=0, 
OF = 0.5 and (2B,t y? = 1.0. The abscissa axis is shown 
E’RC 
as a function of O7/O5 (or T/T Op): The curves asymptote 
O 


to (N-1)/2. Some observations are of interest. First, the 
1 
slope near the origin is equal to (2Bgt pg) OF ee So nCces Lol 


a given quantizer E[I]/a reaches a fixed value for large 
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values of O7/05 » increasing (2Bot 8 (which physically 


2 RC? 
means that more smoothing will be done by the RC filter) 


decreases the range of linear output. Therefore t Should 


RC 
, be chosen as the minimum value which will enable the sampler 


to operate properly. Second, the parameter oF will be 
O 
determined by the amount of degradation tolerated at Top" 


And third,.the desired range of linear output will determine 
the number of steps of the quantizer. Therefore a procedure 


vor use a Figure such as 56 would be as follows: First, est 


. 
at Spm Dot 


1 _ 
Sy and on using the criteria given above. 
O 


second, draw a set of curves of the expected value of the 


determine (2Bot 


output versus temperature using (156). Third, decide upon 
a criterion of how much departure from linearity is allowed 
in the temperature range of operation, and finally obtain 
the required number of steps from the curves drawn. 

To summarize, the effects of the addition of an ADC on 
the linearity of the output of the radiometers under 
consideration, have been investigated in this chapter. The 
knowledge of the behavior of the output is important for 


calibration purposes. 
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values of O7/o5 » increasing (2B,t % (which physically 


x RC? 
means that more smoothing will be done by the RC filter) 


decreases the range of linear output. Therefore t,, should 


RC 
be chosen as the minimum value which will enable the sampler 


to operate properly. Second, the parameter or will be 
O 


determined by the amount of degradation tolerated at Top" 
And third,.the desired range of linear output will determine 
the number of steps of the quantizer. Therefore a procedure 


to use a Figure such as 56 would be as follows: First, 

1 =: 

aU and on using the criteria given above. 
O 


second, draw a set of curves of the expected value of the 


determine (2Bpt 


output versus temperature using (156). Third, decide upon 
a criterion of how much departure from linearity is allowed 
in the temperature range of operation, and finally obtain 
the required number of steps from the curves drawn. 
Lowsummarize, the effects of the addition of an ADC on 
the linearity of the output of the radiometers under 
consideration, have been investigated in this chapter. The 
knowledge of the behavior of the output is important for 


calibration purposes. 


JUS, 





where F is the degradation factor for the equivalent total 
power radiometer. The factor multiplying F is the performance 


figure of a balanced Dicke radiometer. 
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Equation (160) follows from the fact that for the half 
of the switching cycle that the radiometer is connected to 
che source, it behaves as a total power radiometer. 

Care must be taken however on the selection of t. If 
commencement of sampling is delayed after switching, it must 
be corrected for the loss in integration time. Also, for 
bhe case of RC filtering followed by sampling, the introduc- 
tion of the RC filter must be reflected in an increase of T 
since F takes only into account the effects of the inclusion 
of the ADC. 

The problem of dynamic range and linearity of the output 


can be handled with the curves developed in Chapter VI. 
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B. NOISE ADDING RADIOMETER 

Figures 59 and 60 show two digital noise adding radio- 
meters. The one in Figure 59 corresponds to the IF sampling 
case and only one channel is shown. The one in Figure 60 
corresponds to the RC filtering case. The extension of 
results already found to NAR is considerably more difficult 
than for a balanced Dicke radiometer. This is due to the 
processing involved [Ref. 13]. 


Let 


i seE +e (161) 
and 


mL ete (162) 


where i. and I, are the outputs of the radiometer with 


an and Top respectively. En and Bo are the expected 


value of the output for the case of To +t and Top respec- 


p 


tively. En and Ee, are independent, zero mean random variables 


with standard deviation equal to OF and Or and represent 
n O 


the deviation from the mean for a particular measurement 
(one cycle). 


Then the Y-factor is [Ref. 14] 


En 
1, En 1 + = 
Ve=4— =9—"* [ n i) (163) 
IE E ag FE 
Oo O ED . 
1 + 
Oo 


ES 





ef ct and Ee, are small compared with En and Eo? then 


= 78 the - 25 (164) 
oO En Eo 
E[Y] = E (165) 
and z E 
oe = [ i cg? + ge? (166) 


Using Equation (166) and the fact that for a NAR 


th 
7, we 
= (1 + -oP) ( (Le) 
op Th aoe 





where t is the switching period; for the case of Figure 2 


it can be shown by straight substitution that, 








2 
) 
i 7 areI \2 ee “ny 52 , dE[I.] 9 92 
2 O o n 2 ) O 
do O.. do 
xX Ko X 
Eo - ETT J qT, + - ETI] ie 
NAR , a{I_ J 
Oy dE[I J 
(1 + —2)* L doy - — @g2 J 
Oy X 
n LAR | Sk ao me ye 
(168) 


where oO icOewat lame, 0 is )o.2at T 


114 





6S HuNdIA 


ieee 0 cd | 

! ! | ! 

| ee : | - 

: | (WA) 6 C TZ , ! 7 
t ! z 

| =i | | | i 
I 

! | ! (FH —@ 

Ps GAC F Hass Say 

! > I z : = : 

: = x | | 

| | (nb Coz | 

| | ! 

es 


15 





H 


{ 
| 
j 
| 
i 
! 
| 
t 
| 
{ 
{ 
| 
! 
| 


= = —_— = = =a = am ~~w = =) =e a gm == an «= ra Dee 


V=9 





io f 
a 


3 
AN 
gl 





tu 


a 
> 
D aad 
=} 
4 
= 


V 


[ 


WV 


—_— PP Pe ee ee ee 


09 dYNnoId 





3 


ia 


Lf 


— a eee ee eee 


116 





i is the degradation factor of the corresponding total 


power radiometer at TeT tT, » and Bo is the degradation 


p 
factor of the corresponding total power radiometer at Tel 5p° 
For the case of Figure 3, of must be replaced by oo slit [Gabe tc) 
All the terms of (168) can be calculated from equations 
previously developed. 

It must be noted that for the case when the signal is 
filtered before sampling (Figure 60), f' (offset) cannot be 


op 
The problem of dynamic range and linearity is easy to 


zero at T or Bs since the ADC would give zero output. 


handle. From Figures 55 or 56, depending on the system under 
consideration, and (165) a plot can be made relating E[T] 

to temperature. Figure 61 shows how this is done for the 
system of Figure 59. 

In summary, in this chapter the results obtained in the 
six previous chapters for the digital total power radiometers 
Shown in Figures 2 and 3, have been used to determine 
performance of other digital radiometers. The balanced 
Dicke and noise adaiae radiometers were investigated. For 
the balanced Dicke radiometer, there is a simple and straight- 
forward Mel dt ionship. For the NAR, however, a relationship 
was found but it is more involved. The problem of dynamic 


range and linearity can be handled easily. 


ily; 
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In the analysis of the classical total power radiometer 
an integrator is used to smooth the output signal. This 
scheme works well but it has the disadvantage that it gives 
output information only at discrete intervals of time, corre- 
sponding to blocks of data. Such performance might not be 
acceptable in some cases because a smooth output might be 
desirable. Therefore, a pure integrator is not generally 
used but rather a low pass filter (such as a RC) that pro- 
vides a continuous output. The theory and equations of a 
classical total power radiometer do not change when using 
these low pass filters instead of an integrator if an 
equivalent integration time is defined [Ref. ji]. 

The digital total power radiometers of Figures 2 and 3 
suffer the same disadvantage as the classical TPR. That is, 
output information is obtained only after the summers have 
completed the summation of the n samples. Therefore a scheme 
where output information is obtained more often must be 
investigated. This chapter deals with ene substitution of 
the Soe ee ioe Figures 2 and 3 by digital filters. There is 
an additional reason for wanting to do this analysis. Since 
a computer is already being used in the radiometers of 
Figures 2 and 3, the use of a digital filter often does not 
require any more hardware. 

Figure 62 shows a digital filter operating on a digital 


Signal X The filter is characterized by its impulse 


Feu 


le 





FIGURE 62 


response Ay and the output is Tye The filter is assumed 
to be causal (i.e. h,=0 for k<0) and time invariant. This 
filter can replace the summers of the radiometers defined in 
the previous chapters. 

Referring then to Figure 62 


oo 


eee xX, bone (169) 


j=—0 


By straightforward algebra it can be shown that 


Pow, Jes E(Y | stEL XI: Va (170) 
k pe ai 
i1=0 
and 
2 ic = Dg 
Oy = Ny hy 2 [Ry -~ E (X]Jh,. ery 1) 
1=0 m=0 m- | 
where ras = ELK, Xan? is the discrete autocorrelation of X- 
Let 
jzri-m (172) 
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then — 


a al 
ope E£ hy Ff (Ry - E°Cx]]n,_, (173) 
i 8. jz-© 4 
Assume that D, decays much slower than Ry - This can be 


seen by observing the location where this filter would go 
in Figures 2 and 3. For the case of Figure 2, the power 
spectral density at the input of the digital filters is 
shaped by the RF, Mixer, IF filter which is broad-band. 
For the case of Figure 3, it can be assumed only when tac 
is small (this is the case for which Figures 46 through 51 
apply). 

Then (173) becomes 


oo f 
o._ == z h ¥ [Ry 


2 
: gece Xy > E[X]In, (174) 


Let J, be the value where R(j,) is negligible. Then 


(174) can be rewritten as 


1 
ne of (Ry - EEX]! 
er 


dG 1 
z CR 
0 j=-@ 


— acer 2 


im 8 


Ky J 
G75) 


It can be shown that the first summation is much smaller 
than the second, and hence it can be ignored. Furthermore 


the limits of the second summation can be changed as follows 


without altering the value of of significantly 


Ten 





ae hy oe (Ry = E°(x)] (176) 


Mor Lnes Summer Wsed in previous chapters, it can bewshown 


(using the same notation of Figure 62) that 


n-1 n-1 


A 2 2 k D 
Geeta se TR FE (xl |- eas = [Ree =] EO X 1 
and 
Se] eS ea (178) 


The upper limit of the summation of (177) can be replaced 
by infinity (due to the narrowness of R,(u)) and the third 
term can be ignored (see Chapters III and V). 

Substituting (170) and (176) into (4) for the general 


veer eer Tere it follows that 








2 
Oo oO = CO CO 
5 > hs fi 5 rhe [R, ~ Ex] 
= = : -~ aL ‘il 
T i=0 T $=0 j=l *j 
AT.2 © lp 
— =e = 79) 
as | a) del xine T 
( z h, ) ( wie | op 
Doing the same with (177) and (178) it follows that 
Zz 
Oo [2 2) 
X 2 2 
ob free hee ia! - E(x] 
( _, 5 ee eee eee (180) 
op pea)? ths 
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Now if a comparison is made of (179) and (180) it is 


clear that they are equal if n, is defined as follows 


E 


M1 8 


)e 


b> 


lim gi ll 
> 


(181) 


jute 


and n, can be interpreted as an equivalent summation number. 


E 
Physically it is @ measure of the performance of a general 
digital filter in terms of the performance of a pure summer. 
in other words, it provides a way to replace thei digital 


filter by a summer of n_ Samples without changing the per- 


EB 
formance of the radiometer. Then for a given sampling rate 


and known n, an equivalent integration time can be determined 


E 
and a performance figure can be calculated using the curves 
of Chapter V. 

For a digital total power radiometer such as those of 


Figures 2 or 3 using a digital filter, a would be given by 
op 


op ea) 
where F is the degradation factor for the equivalent 
—l, 
radiometer using summers and (Bon,A) * is the performance 
figure of a total power radiometer with integration time 


equal to n,A. 
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It is often convenient to express (181) in its equivalent 
Z-transform since the transfer function of a dae ion erat ex 
is usually given in terms of its z-transform. It can be 
shown that [Ref. 15] 


2 
CHCZ =) 
n. = al (183) 


1 -l, _-l 
On é HUZ JH Zi) 2 de, 
r 
\ 
where the contour integration denoted by [f is around a unit 
cirele and H(z) is the transfer function of the filter. 
. AS an example of the use of (183) a first and second 
Srocreqicitalefiiter will be considered. 


For a first order filter, the transfer function is 


given by [Ref. 13] 


Wea = (184) 


io rz 


It can be easily shown that (183) becomes 





no = I - (185) 


The transfer funétion of a second order filter Is given 


by | 
eo. (186) 


Substituting (186) into (183) and after considerable algebra, 


it can be shown that 
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. 4 a 1© (a~+2) _ 2raé — am + 1 


n (187) 
ze Gee? Gee) 

Equation (182) and (183) are very similar to the ones 
defining the equivalent integration time of a low pass filter 
of an analog total power radiometer. This is not surprising 
Since many similarities exist between digital and analog 


networks even though they are not equivalent. 
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IX. CONCLUSION 


This work investigates the performance of radiometers 
using digital processing. It is assumed that the changing 
of the signal (esac an analog to digital form is performed 
before the final output is obtained. Two digital total 
power radiometers were considered in detail. 

In Chapter II a figure of merit is defined in order to 
evaluate performance. The effects of the addition of an 
analog to digital converter are discussed. 

In Chapter III the effects of an ideal ADC are investi- 
gated. It serves as an introduction to the problem of 
sampling and quantization. The results are useful since 
they provide a lower bound for the degradation factor as 
defined in Chapter II. 

Chapter IV deais exclusively with the effects of a 
quantizer on the autocorrelation function of a zero mean 
Gaussian input process. The problem of a quantizer followed 
by a square law detector is investigated. 

Chapter V combines the information of Chapters II, Iil 
and IV and treats the effects of a real ADC. Both quantiza- 
tion and sampling effects are investigated. For one of the 
radiometers under investigation, exact results are obtained. 
For the other a reasonably good Acorn ments made and an 
approximate solution is obtained. Other approximation are 


discussed and compared. 
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Chapter VI discusses the problem of dynamic range and 
linearity involved in a digital radiometer. 

In €hapter VII the results obtained in the previous 
Chapbers are used to investigate other types of digital 
radiometers. Balanced Dicke and noise adding radiometers 
re eels HaCTCd. 

Chapter VIII investigates the problem of using a general 
digital filter for processing of Signals in a radiometer. 

The major difficulty encountered throughout the work was 
the number of aes Lhae CGivcwe ne emma =p rolb Lem. Ne ree 
fore it was impossible to do a global analysis. However, 
the factors and considerations that enter into the design of 
digital radiometers were treated in sufficient depth so that 
a good understanding of the advantages, problems and 


limitations of digital radiometers is presented. 
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APPENDIX 


The following is a list of the programs used to 


evaluate the equations and draw the figures throughout the 


They are listed in the order they were used. 
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